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I. INTRODUCTION 

This report is concerned with the construction of a high- 
frequency solution for the diffraction of an electromagnetic wave 
obliquely incident on a curved edge in an otherwise smooth, curved, 
perfectly-conducting surface surrounded by an isotropic, homogeneous 
medium. The surface normal is discontinuous at the curved edge, and 
the two surfaces forming the edge may be convex, concave or plane. 

The solution is developed within the context of Keller's geometrical 
theory of diffraction * * (referred to simply as the GTD henceforth) 
so the dlyadic diffraction coefficient is of interest. Particular 
emphasis is placed on finding a compact, accurate form of the diffraction 
coefficient valid in the transition regions adjacent to shadow and 
reflection boundaries and useful in practical applications. 

According to the GTD, a high-frequency electromagnetic wave incident 
on a curved surface with a curved edge gives rise to a reflected wave, 
an edge diffracted wave, and an edge excited wave which propagates 
along a surface ray. Such surface ray fields may also be excited at 
shadow boundaries of the curved surface. The problem is easily 
visualized with the aid of Fig. 1, which shows a plane perpendicular 
to the edge at the point of diffraction Q^. The pertinent rays and 
boundaries are projected onto this plane. To simplify the discussion 
of the reflected field we have assumed that the local interior wedge 
angle is ^ ir. According to Keller's generalized Fermat's principle, 
the ray incident on the edge produces edge diffracted rays ed and 
surface diffracted rays sr. In the case of convex surfaces the 
surface ray sheds a surface diffracted ray sd from each point Q on its path. 




Fig. l. Incident, reflected and diffracted rays and their 

associated shadow and reflection boundaries projected 
- onto the plane normal to the edge at the point of 

diffraction Q^. 

ES is the boundary between the edge diffracted rays and the surface 
diffracted rays; it is tangent to the surface at Q^. SB is the shadow 
boundary of the incident field and RB is the shadow boundary of the 
reflected field, referred to simply as the reflection boundary hence- 
forth. If both surfaces are illiminated, then there is no shadow 
boundary at the edge; instead there are two reflection boundaries for 
the problem considered here. Since the behavior of the ray optics 
field is different in the two regions separated by a boundary, there 
is a transition region adjacent to each boundary within which there 
is a rapid variation of the field between the two regions. 
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4 V In the present analysis it is assumed that the sources and 
field point are sufficiently removed from the surface and the 
boundary ES so that the contributions from the surface ray field can 
be neglected. The total electric field may then be represented as 

(1 j F = e‘* u'' + u’^ + 

in which is the electric field of the source in the absence of 

the surface, F^ is the electric field reflected from the surface with 

the edge ignored, and F^ is the edge diffracted electric field. The 
i r 

functions u and u are unit step functions which are equal to one in 
the regions illuminated by the incident and reflected fields and to 
zero in their shadow regions. The extent of these regions is determined 
by geometrical optics. The step functions are shown explicitly in Eq. (1) 
to emphasize the discontinuity in the incident and reflected fields at 
the shadow and reflection boundaries, respectively. They are not included 
in subsequent equations for reasons of notational economy. 

The diffracted field as defined by Eq. (1) penetrates the shadow 
region, which according to geometrical optics has a zero field, to 
account for the nonvanishing fields known to exist there. But the 
correct high-frequency field must be continuous at the shadow and 
reflection boundaries; hence the diffracted field must compensate 
the discontinuities in the incident and reflected fields there. 

In other words, the diffracted field must provide the correct 
transition between the illuminated regions and the regions shadowed 
by the edge. 
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The high-rfrequency solution described in the next sections is 

4 

obtained in the following way. A Luneberg-Kline expansion for the 
incident field is assumed to be given. The reflected field is 
expanded similarly and related to the incident field by imposing 
the boundary condition at the perfectly-conducting surface. Only the 
leading term, the geometrical optics term, is retained.^,, Next the 
general form of the leading term in the high-frequency solution for 
the edge- diffracted electromagnetic field is determined. The wedge 
(straight edge) geometry is treated first; its dyadic diffraction 
coefficient is deduced from the asymptotic solution of, several canonical 
problems. Some parameters in this diffraction coefficient are seen 
to depend on the type of edge illunination. They are determined for 
an arbitrary incident wavefront by requiring the leading term in , 
the total field to be continuous at the, shadow and reflection boundary. 
It is found that only a slight extension of the solution for. the wedge 
is needed to treat the more general problem posed by the curved edge. 

This report is the third in a series of reports dealing with 

5 

edge diffraction. In the first report , the Pauli-Clemmow method of 
steepest descent was employed in a manner different from that 
employed by Pauli® to obtain a more accurate asymptotic solution for 
the field diffracted by a wedge. We showed that our, generalized Pauli 
expansion can be transformed term by term into a generalized form of 
the asymptotic expansion given by Oberhettinger^. The leading term 
in our expansion was found to be more accurate than the leading term 
in Oberhettinger'sexpansipn; furthermore, our leading term for the 
diffracted field contains a simple correction factor, which permits 
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the field to be calculated easily in the transition region. This 
property is of (bhsiderable practical importance, because it enables 
one to use the geometrical theory of diffraction in the transition 
regions without introducing a supplementary solution. The correction 
factors, referred to here as transition functions are simply included 
with the diffraction coefficient. 

In the first report only the scalar problem of plane waves 
normally incident on the edge of wedge is considered. In the second 
report this work is extended to obtain a dyadic diffraction coefficient 
for a perfectly-conducting wedge illuminated by obi iquely-inci dent 
plane, cortical, and spherical waves. By introducing the natural, 
ray-fixed coordinates, the dyadic diffraction coefficient obtained 
from each of these canonical problems is reduced to the sum of two 
dyads. In other words, the matrix formed by the elements of the 
dyadic diffraction coefficient is a two by two diagonal matrix. 

The diagonal elements of this matrix are simply the scalar dif- 
fraction coefficients, and D^, for the Neumann (hard) and 
Dirichlet (soft) boundary conditions, respectively. The transition 
functions appearing in and have the same form for the four 
types of illumination; in each case only a Fresnel integral is 
involved. However the argument of the Fresnel integral depends 
upon the type of illumination. Outside of the transition regions these 
factors are approximately one, and Keller's expressions for the dif- 
fraction coefficients are obtained. The asymptotic solutiorts 
described in this paragraph help us formulate the solution for a 
more general type of illumination of the wedge, as noted earlier. 
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The analysis of wedge diffraction has had a lengthy history. - 

Only a few of the reports and papers have been mentioned thus far. 

Many of the more important papers on this subject may be found in 

references 9, 10. A good review of wedge diffraction and the special 

case of half plane diffraction is given in Chapters 6 and 8 of 

reference 9. Recently Anluwalia, Boersma, and Lewis have written some 
11 12 13 

papers * * of special relevance to the work described here. References 

11 and 12 describe high frequency asymptotic expansions for scalar waves 
diffracted by curved edges in plane and curved screens and 
reference 13 extends this work to a curved edge in a curved surface. 

The authors make use of ray coordinates, and some of their results 
dealing with rays and wavefronts have been helpful in the development 
of our solution. Nevertheless, there are some noteworthy differences 
between their solutions and this one, apart from the fact that their 
problem is scalar instead of the vector problem treated here. Their 
formulation or ansatz begins with the total field, and the 
resulting correction of the ordinary GTD solution in the transition 
region is different from ours. It appears that our result is the 
more convenient to apply. 

II. THE GEOMETRICAL OPTICS FIELD 

The geometrical optics field, which is the sum of the leading 
terms in the asymptotic expansions for the incident and reflected 
fields, is a part of our high frequency solution for the edge 
diffraction. The asymptotic expansions for the incident and 
reflected fields are presented in this chapter; the results are 
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not new and so they may be familiar to the reader, but they are 
included here for the sake of completeness and continuity in the 
discussion. 

- Let the smooth, curved perfectly-conducting surface S be a 
part of our curved edge structure; this surface is defined by 

the position vector Rg(u,v), where u and v are the curvilinear 

^ A /W A A 

surface coordinates with u • v = u • n = v • n = 0 in which the 

A 

superscript denotes a unit vector and n is the unit vector 
normal to S, see Fig. 2. 



Fig. 2. Geometry associated with the reflection by a 
curved surf ace S. 


A high-frequency electromagnetic wave propagating through an 
isotropic, homogeneous medium is incident on S. The incident and 
reflected fields and quantities associated with them are denoted 
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by the superscripts (or subscripts) i and r, respectively. The 
boundary condition on the electric field at S is 


(2) fixer' (l?j) +r(TTj)] = 0. 


Since we are interested in an asymptotic high frequency solution, 
the incident and reflected fields are expanded in Luneberg-Kline 
series for large u 


(3) 


rof) . I 

m=0 (juj) 


m 


where an e'^“* time dependence is assumed and k = u/Vp with Vp the 
phase velocity of the medium. The electric field is a solution 
of 


(4) (v^ + k^)r = 0 

subject to the condition that 

(5) V . E = 0. 

Substituting Eq. (3) into Eqs. (4) and (5) and equating the 
coefficient of each power of w to zero, one obtains the eikonal 
equation 

(6) |v,J;r=l, 

together with the transport equations 
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whose solutions must also satisfy 
(8a) 4: v^.= .0 ., 

(8b) ^ ^ * ^m-1* ” 1^2,3 ••*, 

S • ’ . 

A 

where v.\p = s a unit vector in the direction of the ray path, which , 
is normal to the wavefront li/CR) = constant, and s is the distance 
along the ray path. 

We are interested here in the solution at the high frequency 
limit, where the asymptotic approximation for T reduces to 

(9) r(s) F^(s) 

Equation (7a) is readily integrated, and after some manipulation 
one obtains^^’^^ 

^o(') = ^o(°)j Tp - , T ])(p,n) 

in which s=0 is taken as a reference point on the ray path and pp 
P 2 are the principal radii of curvature of the wavefront at s=0. 

In Fig. 3 p^ and P 2 are shown in relationship^, to the rays and , 
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wavefronts. Employing the Maxwe ,7 cuH eguatfon v . r = n 

follows from Eq. (3) that the leading tern in th 

« leaaing term in the asymptotic ap- 

proximation for the magnetic field is 


ni) 

c 


where =/^ ts the characteristic 

E is given by Eq. (9). 

Employing £q, (5), 

. s ds = ds; 


admittance of the medium, and 


consequently. 


02 ) 


V'(s) = tp(o) + s. 



Fig, 3. Astigmatic tube of rays 
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Frorh'Eqs. (10), (12) and (9) one obtains the leading term in the 
asymptotic expansion 


(13) 


F(s) 


''' E^(o) 


,-jk'l'(o), 


f P*j P2 
(P'j'^’S ) (p2 




,-jks 


which is ipecognlzed: as the geometrical optics field; this could have 
been deduced from classical geometrical optics employing power 
conservation in the tube of rays shown in Fig. 3. 

It is apparent that when s = -p^ or -p 2 Eq. (13) becomes infinite 
so that it is no longer a valid approximation. The congruence of 
rays at the lines 1-2 and 3-4 of the astigmatic bundles of rays is 
called a caustic. As we pass through a caustic in the direction of 
propagation the sign of p+s changes sign and the correct phase shift 
of +tt/ 2 is introduced naturally. Equation (13) is a valid high 
frequency approximation on either side of the caustic; the field at 

V 

a caustic must be found from separate considerations^^*^^. 

Returning now to our problem of reflection at the. perfectly- 
conducting surface S, the incident field is known thus the 
are known. The found by using the boundary condition. 

We take the: surface Sai the reference point on the reflected ray 
so s = 0 thbre; furthermore let' us modify our notation for, the 
incident field at point oC reflection on S, denoting it by 
r^(Q[^). Now substituting Eq. (3) into Eq. (2) and equating like 
powers of u). 
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(14) n X E^(Qp)e ^ + n x eJJCo)? J = 0, m=04,.2,3...--^y 

If this equation is to be satisfied for all m, then 

(15) iP^(Qr) = H^^(o). 

and since the above is true for every point on S, 

(16a) u • Vi//^.(Qj^) = u • Vt|;^(o), 

(16b) V • V(i/^-(QR) = V • Vij;^(o), 

or ^ : •••57 

(17a) s^ • u = s’^ • u, 

(17b) s^ • V = s*^ • V . 

Thus s^ and s*” have the same projection on the plane tangent to S 
at Qp, which leads us to the law of reflection : 

The plane of incidence is defined by the incident ray and the 
normal to the surface at the point of incidence. The reflected ray 
lies in the plane of incidence and the angle of reflection equals 
the angle of incidence e\ where both angles are measured from the 
normal to the surface as shown in Fig. 2. 
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\ 

\ 

\ — r 

! To (tetermine Eq{o) it is convenient to introduce the unit vector 

\ A 

* ej^ perpendicular to the plane of incidence and the unit vectors e„ 

and ej,, which are parallel to the plane of incidence and perpendicular 

' I 

to s and s , respectively, so that 

(18) X s = e„ 

in each case. Then we may set 


F = e E + e E 
0 II on J- ox 


at the surface. Employing Eqs. (15) and (19) in Eq. (14) we obtain 


Ej(o) = rICQj) . TT= fJ(Qs) • [e>-,r - 


where F is the dyadic reflection coefficient. In matrix notation 
the reflection coefficient has a form familiar for the reflection 
of a plane electromagnetic wave from a plane, perfectly-conducting 
surface 



This is not surprising if one considers the local nature of high- 
frequency reflection, i.e., the phenomena for the most part depends 
on the geometry of the problem in the inmediate neighborhood of Qp. 
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thus the surface S can be approximated by its tangent plane at Q[^, 
and the wavefront of the incident field by a plane wavefront. 

It follows from Eqs. (13), (15) arid (20) that the geometrical 
optics reflected electric field 

(22) rCs) = r'(Q„) • T 

In which p!| and pL are the principal radii of curvature of the 
reflected wavefront at the point of reflection Qj^. In Appendix I 
these radii of curvature are found to be a function of the incident 
wavefront curvature, the aspect of incidence and the curvature of 
S at Qp. In this appendix it is shown that the expressions for 
and pg can be put into the form 



(23a) 



(23b) ^ 

^2 




9 


where and are the principal radii of curvature of the incident 
wavefront at Qj^ and 



A A 

with X-| and X2 unit vectors in the principal directions of the 
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incident wavefront, U-j and, (J 2 unit vectors in the. principal directions 
of S at Qp, 

Equations (23a, b) are reminiscent of the simple mirror formulas of 
elementary physics; this is particularly true ip the case of an 
incident spherical wave, where = s' and f^ f 2 are focal 

distances independent of the range of the source of the spherical 
wave. 

In principle the geometrical optics approximations can be 
improved by finding the higher order terms Elj'(R), T 2 (^» *** in 
the reflected field, but in general it is not easy to obtain these 
from Eqs, (7b), (8b), (14) and (15). Furthermore, these terms do 
not correct the serious errors in the geometrical optics field 
resulting from the discontinuities at reflection and shadow boundaries. 
In the next section we will construct high-frequency approximations 
for the edge diffracted field which in combination with the geometrical 
optics field yield a continuous total field. : 

III. THE EDGE DIFFRACTED FIELD 

The smooth surface S has a curved edge formed by a discontinuity 
in its unit normal vector. Points on the edge are defined by the 
position vector F. When an electromagnetic wave is incident on 
the edge a diffracted wave emanates from the edge. The leading, - 
term in the .high frequency approximation for the electric field is 
assumed to have the form ' •: ' 

(24) A(^) . 


and R-|,R 2 the principal radii of curvature of. S at Qp. 
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Substituting the above expression for into Eqs. (4) and (5), and 
again equating the coefficient of each power of co to zero, one obtains 


(25) 


= 1 . 




^ A + 1 (S^rp)A = 0, 


(27) 


s-A = 0, 


and from the discussion in the preceding section, it follows that 


(28) 

in which 


E^ 




( 29 ) . 

JT 

and s is the distance along the diffracted ray from a reference point 
O' which is not a caustic of the diffracted ray, see Fig. 4. 

It is convenient to locate the reference point of the diffracted 
ray at the edge point from which it emanates; however the edge is 
a caustic of the diffracted field. On the other' hand, it is clear 
that ^(s) given by Eq. (28) must be independent of the location 
of O' , hence 

lim 1^(0' )V^ exists. 

Furthermore E^(s) is proportional to the incident electric field at 
p£, so we may set 
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(30) lim E^Qp) • D , 

p'.-^ : ■ “^ • 

where ^ is the dyadic edge diffraction coefficient, which is analogous 
to the dyadic reflection coefficient of the preceding section. 



Thus the edge diffracted electric field 


in which p is the distance between the caustic at the edge and the 
second caustic of the diffracted ray. 

In appendix II it is shown that 


( 32 ) 



rig • (s'-s) 
a sin^e^ 
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wherein is the radius of curvature of the incident wavefront at 

/s 

taken in the plane containing the incident ray and e is the unit 

A ' y 

vector tangent to the edge at Q^, n^ is the associated unit normal 
vector to the edge directed away from the center of curvature, a>0 
is the radius of curvature of the edge at Q^, and 3- is the angle > 
between the incident ray and the tangent to the edge as shown in 
Fig. 5a. Equation (32) is seen to have the form of the elementary 
mirror and lens formulas in which f is the focal distance. If p is' 
positive, there is no caustic along the diffracted ray path; however 
the caustic distance p is negative if the (second) caustic lies 
between and the observation point. The diffracted field calculated 
from Eq. (31) is not valid at a caustic, but as one moves outward from 
Qg along the diffracted ray, a phase shift of +ir/2 is introduced 
' naturally after the caustic is passed as in the case of; the geometrical 
optics field. 

Since the high frequency diffracted field has a caustic at the 
edge Eq. (31) is not valid there, and we can not impose a condition at 
Qg to.determine "D in a manner similar to that used .to find r". Never- 
theless, the matching of the phase functions at the edge 

(33) 

is a necessary condition, which yields some useful information about 
the solution. Since the phase is matched at each point on the edge, 
it follows from Eq. (33) that 

e • Vi|;.(Q^) = e • Vi|;^(Q^) :=;e r Vi);j(Q^.), 
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i.e., 

^ A A|f% A A 

(34):' ' e • s: = e*s =e*s. 

The angle of incidence in this case is 3^ defined earlier and shown 
in Fig. 5a, cos Bq = e • s^ , 0 I 3^ ^ tt/ 2. The angle of diffraction 
is the angle between the diffracted ray and the tangent to the 
edge at Q^; cos 3^ = e • S, 0 ^ 1. Keller's law of edge dif- 

fraction follows from Eq. (34). 

The law of edge diffraction ; the angle of diffraction 3^j is equal 
to the angle of incidence 3^, so that the diffracted rays emanating 
fromiQg form a. cone whose half angle is 3^ and whose axis is the 
tangent: to the edge. The incident ray and the ray reflected from 
the 'Surface at also lie on the cone of the diffracted rays. 

The form of the dyadic diffraction coefficient will be treated 
next. If an edge-fixed coordinate system is used to describe the 
components of the incident and diffracted fields, it has been found 
that'the dyadic diffraction coefficient is the sum of seven dyads[18,19]; 
in matrix form this means that the diffraction coefficient is a 3 x 3 
matrix with 7 non-vanishing elements. However from Eqs. (8a) and (27) 
it is apparent that if a ray-fixed coordinate system were used in place 
of the edge-fixed coordinate system, the diffraction coefficient would 
reduce to a 2 x 2 matrix, so that no more than four dyads would be 
required. A further reduction in the number of dyads can be anticipated 
if the proper ray- fixed coordinate is chosen. Recall that this kind of 
simplification is achieved in the case of the dyadic reflection coefficient. 
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A 



(b) 


Fig. 5. 


20 



if the incident and reflected fields are resolved into components 
parallel and perpendicular to the planes of incidence and reflection, 
respectively, where the plane of reflection, which contains the normal 
to the surface and the reflected ray, coincides with the plane of 
incidence. Analogous planes of incidence and diffraction can be defined 
in the;; present case; 

The plane of incidence for edge diffraction, referred to simply 
as the edge-fixed plane of incidence henceforth, contains the incident 

A 

ray and the unit vector e tangent to the edge at the point of incidence 
Q^. The plane of diffraction contains the diffracted ray and e. 

These planes are depcited in Fig. 5 ; they are azimuthal planes with 
respect to the polar axis containing e, and their positions can be 
specified by the angles and (j> shown in Fig. 5 b. The unit vectors 
and I are perpendicular to the edge fixed plane of incidence and 
the plane of diffraction, respectively. The unit vector s = s' is 
in the direction of incidence at the edge and the unit vector s is in 
the direction of diffraction. The unit vectors 3 V and 3^ are parallel 
to the edge fixed plane of incidence and the plane of diffraction, 
respectively, and . r 

( 35 a, b) 3 q = s' X , Sq = s x 

Thus the coordinates of the diffracted ray (s,ir-3Q,4)) 
spherical coordinates and so are the coordinates of the incident ray 
(s' ,3 q,())' ), except that the incident (radial) unit vector points toward 
the origin Q^. 
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According to Keller's theory[3] the diffraction coefficient for 
a curved edge may be deduced from a two-dimensional canonical problem 
involving a straight edge, where the cylindrical surfaces which form 
the edge are defined by the boundary curves depicted in Fig. 5b, In 
the present discussion the edge may be an ordinary edge formed by a 
discontinuity in the unit normal vector, an edge formed by a discontinuity 
in surface curvature, or an edge formed by a discontinuity in some; higher 
order derivative of the surface. 

Consider the z-components of the electric and magnetic fields; in; 
the presence of this surface with an edge 

(36a) ^ 


(36b) 


they satisfy 


(37) (v^ 




"z 

^ J 


r= 0 


together with the soft (Dirichlet) or hard (Neumann) boundary 
conditions 


(38.39) Ez = 0 O'" = 0. 

respectively, on the boundary curve and the radiation condition 
at infinity. The a/ an is the derivative along the normal to the 
boundary curve. 
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starting with the high frequency solutions for the z-components 
of the diffracted field 


g- ■ ^ A(R), 


and ;substituti rig these into Eq. (37), and employing the methods 
described earlier, the asymptotic solutions may be put into the form 


4 4 


4 “h 


S(p+S, 


in which is referred to as the soft scalar diffraction coefficient 
s 

obtained when the soft boundary condition is used, and is referred 
to as the hard scalar diffraction coefficient obtained when the 
hard boundary condition is used. 

Since 


(42a) e1 = e\ sin 6„ . 


(42b) e;, sin 




(43a) E^ = -E^ sin e^. 


(43b) ■ -V^ E^ sin s„. 


with 1/Y^ = =y p/e the characteristic impedance of the medium. 
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it follows from Eqs. (41), (42), and (43) that 


(44) 


>= - 


“h 


N 


P .-Oks . 

Hp+iT 


consequently, the dyadic diffraction coefficient for an ordinary (or 
higher order) edge is a perfectly-conducting surface can be expressed 
simply as the sum of two dyads 


(45) D = 

to first order. Since Dg and Dh are the ordinary scalar diffraction 
coefficients which occur in the diffraction of acoustic waves which en- 
counter soft or hard boundaries, we see the close connection between 
electromagnetics and acoustics at high frequencies. Also, it follows that 
the high frequency diffractiop by more general edge structures, and by thin 
curved wires can be described in the form given by Eqs. (44) and (45). 

The balance of this report is concerned with finding expressions 
for Dg and which can be used in the transition regions adjacent to 
shadow and reflection boundaries in the case of diffraction by an 
ordinary edge. Recently Keller and Kaminetzky[20] and Senior[21] 
have obtained expressions for the scalar diffraction coefficients in 
the case of diffraction by an edge formed by. a discontinuity in 
surface curvature and Senior[22] has given the dyadic (or matrix) 
diffraction coefficient in an edge- fixed coordinate system. Keller 
and Kaminetzkey[20] also have given expressions for the scalar 
diffraction coefficients in the case of higher order edges. 

The diffraction by a wedge will be considered first; the straight 
edge serves as a good introduction to the more difficult subject of 
diffraction by a curved edge. As noted earlier, the dyadic diffraction 
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coefficient can be found from the asyn^jtotic solution of several 
canonical problems, which involve the illumination of the edge by 
different wavefronts. It is not difficult to generalize the resulting 
expressions for the scalar diffraction coefficients to the case of 
illumination by an arbitrary wavefront. 

A. The Wedge 

When a plane, cylindrical, or conical electromagnetic wave is 
incident on a perfectly-conducting wedge, the solution may be 
foTTnuTated in terms of the components of the electric and magnetic 
field parallel -to the edge; we will take these to be the z-components. 

In the case of a spherical wave it is convenient to use the z-components 
of the electric and magnetic vector potentials. These z-components 
may be -represented by eigenfunction series obtained by the method of 
Green's functions. The Bessel and Hankel functions in the eigenfunction 
series are replaced by their integral representations and the series' 
are then summed leaving the integral representations. Integral 
representations for the other field components in the edge-fixed 
coordinate system are then found from the z (or edge)-components , 
except in the case of the incident spherical wave, where the integral 
representations of the field components are obtained from the z-components 
of the vector potentials. These integrals are approximated asymptotically 
by the Pauli - Cl eiimow method of steepest descent[23], and the leading 
terms are retained. The field components are then transformed to the 
? ray-fixed coordinate System described previously. The resulting expression 
for the diffracted field Can be written in the form of Eq. (31) which 
makes it possible to deduce the dyadic diffraction coefficient D^. 
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The asymptotic solutions outlined in this paragraph are presented in 
detail in reference [8], 

Summarizing the results given in. Reference [8] 

(46) E^(s) 'v. • f(s\s' ) A(s)e"J''^ 

in which A(s) describes how the amplitude of the field varies along 
the diffracted ray; 

^ for plane, cylindrical and conical wave 
incidence (in the case of cylindrical wave 
incidence, s is replaced by r = s sin 

(47) A(s) =1 the perpendicular distance to the edge), 

for spherical wave incidence. 

It follows from Eq. (32) that p = Pg for the wedge. In the case of 
plane, cylindrical and conical waves is infinite and in the case 
'' of spherical waves = s'. The dyadic diffraction coefficient 

'D(s,s') has the form given in Eq. (45), which supports the assumptions 
leading to that equation. 

if the field point is not close to a shadow or reflection 
boundary, the scalar diffraction coefficients 
(49) ' 

' h ' ^ 

for all four types of illumination, which is important because the 
diffraction coefficient should be independent of the edge illumination 
away from shadow and reflection boundaries where the plane surfaces 
forming the wedge are ({>=0 and cj) = nir. The wedge angle is (2-n)ir; see 
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Fig. 5b. This expression becomes singular as shadow or reflection 
boundaries are approached, which further aggravates the difficulties 
at these boundaries resulting from the discontinuities in the incident 
or reflected fields. The above scalar diffraction coefficients have 
been given by Keller[3]. 

Grazing incidence, where (|)' = 0 or nir must be considered separately. 
In this case = 0, and the expression for given by Eq. (4) must be 
multiplied by a factor of 1/2. The need for the factor of 1/2 may be 
seen by considering grazing incidence to be the limit of oblique 
incidence. At grazing incidence the incident and reflected fields merge, 
so that one half the total field propagating along the face of the wedge 
toward the edge is the incident field and the other half is the reflected 
field. Nevertheless in this case it is clearly more convenient to regard 
the total field as the "incident" field. The factor of 1/2 is also 
apparent if the analysis is carried out with <))' = 0 or nir. 

Combining Eqs. (31), (45) and (49) it is seen that the diffracted 
field is of order with respect to the incident and reflected fields. 

At high frequencies this means that the diffracted field is in general 
weaker than the incident and reflected fields, at aspects not close to 
shadow and reflection boundaries. 

To simplify the discussion, the wedge angle has been restricted so 
that 1 < n ^2; however, the solution for the diffracted field may be 
applied to an interior wedge where 0 < n < 1. The diffraction coef- 
ficient vanishes when 

sin '^ = 0; — 
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hence for n = 1, the entire plane, 

n = 1/2, the interior right angle, 
n = 1/M, M = 3,4,5 ••*, interior acute angles, 
the boundary value problem can be solved exactly in terms of the in- 
cident field and a finite number of reflected fields, which may be 
determined from image theory. Moreover as n -► 0, even with the 
presence of a non-vanishing diffracted field, the phenomenon is 
increasingly dominated by the incident and reflected fields. . . 

Returning now to the subject of exterior edge diffraction, the 
total field changes rapidly in the vicinity of shadow and reflection 
boundaries. In the case of the shadow boundary its behavior is 
predominantly that of the incident field on the illianinated side, 
whereas it is that of the diffracted field, emanating: from' the:edge, 
on the shadow side. For example if the wedge is illuminated by a plane 
wave perpendicular to its edge, the total field varies from an 
essentially plane wave behavior to a cylindrical wave behavior in 
the vicinity of the shadow boundary. These regions of rapid field 
change adjacent to the shadow and reflection boundaries are 
referred to as transition regions . In the transition/regions the 
magnitude of the diffracted field is comparable with the incident 
or reflected field, and since these fields are discontinuous at 
their boundaries, the diffracted fields must be discontinuoTJS at 
shadow and reflection boundaries for the total field to be continuous 
there. 
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in which N" are the integers which most nearly 
satisfy the equations 

(53a) 2iTnN -((j)±({)') = ir 

and 

(53b) ZirnN -((j)±(t>') = — rr . 

The above expression for the soft (s) and hard (h) diffraction 
coefficients contains a transition function F(X) defined by Eq. (51), 
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where it is seen that F(X) involves a Fresnel integral. The 
magnitude and phase of F(X) are shown in Fig. 6, where X = kLa. 

(54a) |F(X)|^1. 

(54b) 0 £ phase F(X) £ tt/4 . 

When 


(55) X > 10, 

F(X) 2^1 . , ■ 

If the arguments of the four transition functions in Eq. (51) exceed 
10, so the transition functions can be replaced by unity, Eq. (50) | 

f 

reduces to Eq. (49). 

(56) X = kL a-(f±<|)' ) 


in which L is a distance parameter, which was determined for several 


types of illumination. It was found that 

2 

s sin for plane wave incidence, 
,1 0 


(57) 


L = < 


r r 


r + r 
ss' 


r for cylindrical wave incidence. 


.•_ 2 , 


5 + jr sin 6^ for conical and spherical wave 


incidence, 

where the cylindrical wave of radius r' is normally incident on the 


edge, and r is the perpendicular distance of the field point from 
the edge. A more general expression for L, valid for an arbitrary 


wavefront incident on the straight edge, will be determined later. 
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The largeness parameter in the asynptotic approximation used to 
find is kL. For incident plane waves the approximation has 
been found to be accurate if kL > 1.0, unless n is close to one, 
then kL should be > 3. 

a“((j) ± 4 )') is a measure of the angular separation between the 
field point and a shadow or reflection boundary. The + and - super- 
scripts are associated with the integers and N“, respectively, 
which are defined by Eqs. (53a, b). For exterior edge diffraction 
= 0 or 1 and N" = -1, 0 or 1. The values of N~ as functions of 
n and 6 = <|) ± 4 )* are depicted in Figs. 7a and 7b; these integers are 
particularly important near the shadow and reflection boundaries 
shown as dotted lines in the figures. It is seen that N" do not 
change abruptly with aspect 4 > near these boundaries, which is a 
desirable property'.' The trapezoidal regions bounded by the solid 
straight lines represents the permissible values of b for 0 £ 4 >, 

4 )' _< nir with 1 1 n _< 2. 

At a shadow or reflection ’boundary one of the cotangent functions 
in the expression for giver, by Eq. (50) becomes singular; the other 
three remain bounded. Even though the cotangent becomes singular, its 
product with the transition function will be shown to be bounded. 
However let us first note the location of the boundary at which each 
cotangent becomes singular; this information is presented compactly 
in Table 1. The locations of the shadow and reflection boundaries when 
only one surface of the wedge is illuminated and when both surfaces of 
the wedge are illuminated are shown in Figs. 8a, b, c below Table 1. 
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TABLE 1 



The cotangent is singular 
when 

value of N 
at the boundary 


4>~. <{* ' ~ n » a SB 
surface 4>=0 is shadowed 

= 0 


(()= (j) ’ + TT , a SB 
surface <{)=mr is shadowed 

N" = 0 

BBS 

(|) = (2n-l , a RB 
reflection from surface (j)=nTT 

= 1 

BiBI 

(p=ir - <j) ' , a RB 
reflection from surface (|)=0 

N" = 0 




Fig. 8. Shadow and reflection boundaries for different 
angles of incidence ((>'. 
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Since discontinuity in the geometrical optics field at a shadow or 
reflection boundary is compensated separately- by one of the four terms 
in the diffraction coefficient, there is no problem in calculating the 
field when two boundaries are close to each other or in juxtaposition. 
This occurs when ({)• = 0 or rnr and when (|)' is close tomr/2 with n 
The shadow and reflection boundaries are real if they occur in physical 
space, which is in the angular range from 0 to mr; outside this range 
they are virtual boundaries. If a virtual boundary is close to the 
surface of the wedge, as it is when (}>' is close to it or (n-l)Tr, its 
transition region may extend into physical space near the wedge and 
significantly effect the calculation of the field there. The value 
of or N“ at each boundary is included in Table 1 for convenience; 
as noted earlier, this is a stable quantity in the transition regions. 

Next it will be shown that the product of the cotangents and the 
transition functions in Eq. (50) is finite, even at the shadow and 
reflection boundaries. To facilitate the discussion let 

(58) 0 = (fi ± (|)' . - " 

In the neighborhood of the shadow or reflection boundary 

(59) 3 = 2tt nN" + (n - e), 

where e is positive in the region illuminated by the incident or 
reflected field. The ± superscript of N is directly associated with 
the + sign in the equation above and the ± sign in the argument of 
the cotangent below. 
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For e small , 

(eojcotfijp). la. , ; 

and 

(61) a-(6) = 2 

The transition function F(X) is given by Eq. (51) with X = kLa~(3).: 

We are concerned with the case where kL is large but X is small, so that 


■‘S) .■(»”) . 

y 


From Eqs. (60), (61) and (62), 


(63) cot^ F[kL a-(6)] 

I f— J H 

nl j2i:kL sgn e - 2kLee je 

for e small. It is clear that the above expression is finite but 
discontinuous at the shadow and reflection boundaries. These 
discontinuities compensate the discontinuity in the incident or 
reflected field at these boundaries, as will be shown in the 
paragraphs to follow. 

The high frequency approximation for the total field being con- 
sidered here is the sum of the geometrical optics field and the 
asymptotic approximation of the diffracted field. It is convenient 
to give the components of these fields in the ray fixed coordinate 
system described on page 21 ; hence it will be necessary to transform 


36 



the components of the reflected field given in the first section to 
this coordinate system. We will begin by carrying out this trans- 
formation, which is facilitated by employing matrix notation. 

From Eqs. (21) and (22) the reflected electric field 



where the subscripts ii and i denote components parallel and per- 
pendicular to the ordinary plane of incidence and 



Note that for the plane surfaces forming the wedge, pjj* = pp = P 2 * 
where p], p^ are the principal radii of curvature of the incident 
wavefront at the point of reflection. Equation (64) may be written 
more compactly as 

(66) e’^ -V re'' f(s) . 

The ordinary plane of incidence and the edge- fixed plane of incidence 
intersect along the incident ray passing through Q^. The ordinary 
plane of incidence, the edge-fixed plane of reflection, and the cone 
of diffracted rays intersect at the ray reflected from Q^. The edge- 
fixed plane of reflection contains the tangent to the edge and the ray 
reflected from Q^. These planes and their lines of intersection are 
depicted in Fig. 9. 
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Fig. 9. Edge fixed plane of incidence and reflection. 

Let the angle between the edge-fixed plane of incidence and the 
ordinary plane of incidence be -a. In Appendix III it is shown that 
the angle between the edge-fixed plane of reflection and the ordinary 
plane of incidence is a. The components of the incident electric field 
parallel and perpendicular to the edge-fixed plane of incidence are 


(67a) e],, = E^ cos a - eJ sin a 

(67b) e\ •= E.’’, sin a + e] cos a 

or in the more compact matrix notation 
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(68) E^' = T(-a)E^' 


where 


(69) T(-a) = 


COSa -Sind 
sincx COSa 


From Eq. (66) the reflected electric field 


(70) e’" -x- R eV/(s) H(e) 

in the neighborhood of the reflection boundary 


(71) H(e) = } (1 + sgn e) 
is the unit step function. 

The components of the reflected field parallel and perpendicular to 
the edge-fixed plane of reflection are given by 

(72) T(a)E^ = [T(a)RT(-a)’''][T(-a)E^'] f(s) H(e) 

From Eq. (69) and R as given in Eq. (21), 

(73) T(a)R T(-a)"^ = R; 

hence from Eqs. (68), (70), (72) and (73) 
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The diffracted field close to the reflection boundary at 4> = Tr- 4 >' 
is given by Eq. (31) together with Eqs. (501 and (63) 



+ terms which are continuous at this boundary. 


For the total field to be continuous at the reflection boundary, 
the sum of the discontinuous terms in Eqs. (74) and (75) must vanish; 
hence 

(76) - 

so that the distance parameter 

(77) L 


s(p^ + s) pj P 2 sin^e^ 

~~ r " i« I I ■ , I I ■ I .1 -I I 

pg(p] + s)(p 2 + s) 



The behavior of the incident and diffracted fields at the shadow 
boundary (|) = tt + (j)' may be treated in the same manner. . After passing 
beyond Q^, the electric field of the incident ray in the neighborhood 
of the shadow boundary is 



The diffracted field close to this shadow boundary is 
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(79) 


ppcn 

^0 


Le 


<P 


1 

'V 2 


1 0 
LO IJ 


rpi 1 


ST 


LeIJ ^oJ s(pI + s) 


sgn e + 


+ in terms which are continuous at this shadow boundary. 


For the total field to be continuous at the shadow boundary, the sum 
of the discontinuous terms in Eqs. (78) and (79) must vanish, and again 
it is' seen that L is given by Eq. (77). Equation (77) is also obtained 
when the leading term in the high frequency approximation for the total 
field is made to be continuous at the other shadow and reflection 
boundaries. Also Eq. (77) reduces to Eq. (57) for the several types 
of incident waves for which formal asynptotic solutions were derived. 

We conclude therefore that the expression for L given by Eq. (77) is 
correct when the wedge is illuminated by an incident field with an 
arbitrary wavefront whose principal radii of curvature are and p^. 

Since kL is the large parameter in the asymptotic approximation, 

3 ^ can not be arbitrarily small, which precludes grazing and near grazing 
incidence along the edge. 

The commentary on Eq. (49) in the case of grazing incidence along the 
surface of the wedge also applies to Eq. (50), i.e., the diffraction 
coefficient is multiplied by a factor of 1/2 and the diffraction 
coefficient = 0. 

If n = 1 or 2, it is apparent from Eq. (52) and the integral values 
of N" that 

(80) a-(3) = a(s) = 2 cos^ . 

Thus 
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(81) cot F[kLa'^(e)] + cot F[kLa"(e)] 

-2 sin ^ 

= = \ F [kLa(3)] , 

cos ^ - cos ^ 

and from Eq. (50) the expressions for the scalar diffraction coefficients 
reduce to 


(82) Dg (<|..<|.';eo) = 

h 


• JL 

4 . ir 

^ n 

nj^ sin 3- 


F[kLa(d)-<f.' )] 

cos ^ - cos 
n n 


F[kLa ((}.+(})' )] 


ir 

COS — 

n 


cos 


•) 


The edge vanishes for n = 1 and the boundary surface is simply a perfectly- 
conducting plane of infinite extent. It is seen that the diffraction 
coefficients and diffracted field vanish for this case as expected. 

If n = 2 the wedge becomes a half plane and \ 


(83) Dg((|),(|)';SQ) = 
h 



Zy/Zirk sin 3 q 


F[kLa(.t-f ) 1 ^ 
cos 


F[kLa(j)+(j)')] 

cos 


which can be written in the form 
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(84) ;3jj) = 

h 


i - 

J /I 


-e 

sT 


'' 4 rrf j2kL cosf 

e V ysgn(TT+ 4 )'- 

j2kL cos^ ^ "I 

sgn(TT-(i)' -(|)) , 


:f f(kL,(j)+<t)' )e 


where 


(85) f(kL,e) = 


.-JT 


dx 


72kL|cos |-| 


» 


a Fresnel integral. 

When the diffraction coefficients given by Eq. (84) are used to 

calculate the fields diffracted by hard or soft half planes il- 

2 

luminated by a plane wave, L = s sin 6^ and the result is in agreement 
with a solution obtained by Sommerfeld[9,24]. Since Sommerfeld's 
solution is an exact solution, we know that our solution is exact 
for this case too. If these half planes are illuminated by a 
cylindrical wave whose radius of curvature is r' , L = rr'/(r+r') in 
which r is the perpendicular distance from the field point to the 
edge, and our solution reduces to an approximate solution deduced by 
Rudduck[25] from the work of 0bha[26] and Nomura[27]. Rudduck and his 
coworkers have applied this solution to a number of two-dimensional 
antenna and scattering problems with good accuracy. 
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In this section on diffraction by wedges, diffraction coefficients 
have been obtained which may be used at all aspects surrounding the 
wedge, including the transition regions adjacent to shadow and 
reflection boundaries. The diffraction by curved edges in plane 
surfaces and curved sheets will be considered in the next section. 

B. The Curved Edge 

The diffraction by curved edges will be treated in this section. 

As in the preceding section our solution is based on Keller's method 
of the canonical problem. The justification of the method is that 
high-frequency diffraction like high-frequency reflection is a local 
phenomenon, and locally one can approximate the curved edge geometry 
by a wedge, where the straight edge of the wedge is tangent to the 
curved edge at the point of incidence in Figs. 5a, b, and its plane 
surfaces are tangent to the surfaces forming the curved edge. The 
reflection coefficient for the curved surface derived in Section I 
could have been found by this method, choosing the reflection of plane 
waves at a plane surface as the canonical problem. With these 
assumptions, the results of the preceding section can be applied 
directly to the curved edge problem. As we have just noted, there is 
an equivalent wedge (with exterior wedge angle nir) associated with every 
curved edge structure, and so in generalizing the solution for the wedge, 
it is only necessary to modify the expressions for the distance parameter 
L, which appear in the arguments of the transition functions. 

In the present treatment we do not show that our solution can be 
matched to a boundary layer solution valid at and near the curved edge. 
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It would be desirable to carry this out to confirm the validity of our 

solution and possibly to obtain additional terms in the asymptotic 

1 3 

approximation. Ahluwalia has used a boundary layer solution in this 
way to obtain an asymptotic expansion for the scalar field diffracted by 
a curved edge;, however his representation of the total field differs 
from the one given here. It does not appear as separate contributions 
from the incident, reflected and diffracted fields. 

The diffraction by a curved edge in a plane screen affords the 
simplest example of curved edge diffraction. The scalar diffraction 
coefficients appearing in Eq. (45) are given by Eqs. (83) or (84), and 
since p = on both the shadow and reflection boundaries, L is the 
distance parameter given by Eq. (77). At aspects other than incidence 
and reflection, p within the square root term of Eq. (31) must be 
found from Eq. (32). As in the case of the wedge, we obtain a high- 
frequency approximation at all points surrounding the edge, which 
are not too close to the edge or to caustics of the diffracted 
field. 

The diffraction by a curved edge in a curved screen (n=2) is 
next in the order of increasing difficulty. Whenever the surface 
forming the edge is curved, the region near it is dominated by surface 
diffraction phenomena, which is particularly important on the convex 
side. On the convex side of the curved screen there are surface 
ray modes, also known as creeping waves, which shed energy tan- 
gentially as they propagate along the surface. As a result of this, 
the radiation leakage phenomenon is significant in a considerable region 
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near the surface. On the concave side of the curved screen we have 
bound modes that do not leak energy as they propagate; these modes' 
are known as whispering gallery modes. Both types of modes are 
excited by an illuminated edge in a curved surface; however they also 
may be excited by the incident field. As mentioned earlier, surface 
diffraction phenomena have been neglected in the present treatment; 
hence the region between the convex surface and the boundary ES 
between the edge diffracted and surface diffracted rays must be ex- 
cluded. The boundary ES is formed by the intersection of the cone of 
diffracted rays and the plane tangent to the surface at Q^; in general 
it does not lie in the ordinary plane of incidence. In addition, the 
transition region adjacent to the boundary must be excluded. This 
region from which the field and source points are to be exluded appears 
as the shaded portion of Fig. 10a, where all rays and boundaries are 
shown projected on the plane perpendicular to the edge at Q^. It should 
be noted that in general the projection of the surface ray does not 
coincide with the intersection of the boundary surface S and the plane 
of projection. 

On the concave side the whispering gallery effect can be 
described approximately by geometrical optics in the form of a series 
of reflected waves whose rays form cords along the concave reflecting 
surface as indicated in Fig. 10b. As glancing incidence is approached, 
the cord length diminishes and the description of the phenomenon in 
terms of a sequence of reflections breaks down; the geometrical optics 
analysis must be truncated at this point. If the errors resulting 
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from this truncation are not serious, the radiation from the concave 
side. can be included in the present analysis. 


A 



\ sr 

SB 


(a) CONVEX SIDE 



(b) CONCAVE SIDE 


Fig. 10. Diffraction at the edge of a curved screen. 

In this case n = 2, and the scalar diffraction coefficients in Eq. 
(45) are given by 
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in which the first term is discontinuous at the shadow boundary, whereas 
the second is discontinuous at the reflection boundary. Unlike the 
reflection from a plane surface, the divergence or spreading of the wave 
reflected from a curved surface is different from that of the incidence 
wave; hence the radii of curvature of the reflected and diffracted wave- 
fronts at the shadow boundary are distinct from the radii of curvature 
of the incident and diffracted wavefronts at the shadow boundary. 
Employing arguments similar to those used to find the distance parameter 
for the wedge 


(87a) 


^ s(pg+s) p]p2 sin^g^ 
Pe(pl'''sXp2'''5) 


, , r s(p^+s) pfpj sin^g 

(87b) l’' = — y-- 2- 

Pe(pi'''s)(p2+s) 


where Pg, pj, p^ are defined as before, plj* and pg are the principal 
radii of curvature of the reflected wavefront at Q^, and from Eq. (32) 


( 88 ) 


A A A A 

1 _ 1 , 2(n.ng)(s'«n) 

T T Z 

*^e ^e a sin'^gg 
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incidence as shown in 


Ac* 4.® a c'Ucic* i.m ar\r\v«/\a 

rw y M mc uppiwuv^ii ^la^iii^ 

M M M 

Fig. 11. Then since q ^ q 2 0> L ->0 and Eq. (86) can no longer be 
used to calculate the scalar diffraction coefficients. Under these 
circumstances the shadow and reflection boundaries usually lie within 
the shaded region in Fig. 11, and the transition regions associated 


with edge diffraction overlap those associated with surface diffraction. 
If the field and source points are both sufficiently far from the edge, 
we may set the transition functions in Eq. (86) equal to unity. On 

. i. 

the otherhand, for the field point or source point close to the edge 
or for both points close to the edge, we may be able to use reciprocity 
(see Appendix IV) to calculate the field at P in Fig. 11, if the distance 
parameters for a unit source located at P are large enough. 



— SHADOW AND 


REFLECTION BOUNDARY'S 
FOR A SOURCE AT P 


Fig. 11. Grazing incidence on the edge of a curved screen. 
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We conclude this chapter by finding the scalar diffraction coef- 
ficients for a curved edge in an otherwise smooth curved surface. Again 
we seek diffraction coefficients which can be used in the transition 
regions associated with the shadow and reflection boundaries of this 
structure. Both surfaces forming the curved edge may be convex, both 
surfaces may be concave, one surface may be convex and the other concave, 
or one surface may be plane and the other convex or concave. 

First let us consider the simple case which occurs when the 
illuminated surface forming the curved edge is plane, as it may be at 
the base of a cylinder or cone. For this configuration the reflected 
field is found directly from the incident field, as it is in the case 
of the wedge, e.g., it may be easily deduced from image theory. Thus 
the scalar diffraction coefficients are found directly from Eq. (50) and 
the distance parameter from Eq. (77). The calculated diffracted field may 
not be accurate close to the shadowed surface if surface diffraction 
phenomena are significant. 

The more general problem where the illuminated surface is curved 
is closely related to the diffraction by a curved edge in a curved 
screen which has just been discussed; for example, the field point 
and source point must not be too close to a convex surface and the 
case of grazing incidence must be treated separately. 

We introduce the wedge tangent to the boundary surfaces of the 
curved edge at Q^. The boundary ES is formed by the intersection of 
this wedge with the cone of diffracted rays. Away from the boundary 
ES on the cone of diffracted rays the scalar diffraction coefficients 
are given by Eq. (50), except that distance parameter L in the argument of 
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each of the four transition functions may be different. As before L 
is found in each case by requiring the total field to be continuous 
at each shadow and reflection boundary. 

It is seen from Figs. 7a, b that N^» N“ associated with the shadow 
boundaries at (|i'-Tr, (j>'+7r are different from zero only at angular 
distances greater than tt from these boundaries. When this angular 
distance exceeds tt the field point is usually outside the transition 
region in question, unless kL is small. In view of the assumptions 
involved in extending the wedge solution to the curved edge, the validity 
of the approximation is in question for such small values of kL, so 
they are excluded here. These considerations and analogous con- 


siderations lead us to set the N“ equal to the values they have in 
Table I. 

Then . v 

. . 'I 

X 


(89) D ) =-^y= 

h ° ZnJZirk sin 6^ 

r 


2 sin F [kL^a(<}.-<|)')] 


cos — - cos 
n 


cot 


F[kL"'a* 

2 8 


( 4 ,+(j)' )]+ cot 




FtkL™a(*+V)] 


in which a(g) = 2 cos 


and a^(3) = 2 cos^ 2rrn^ 

Again employing arguments similar to those used to find the distance 
parameters for the edge, one finds that is given by Eq. (87a), and 
that L*^, are given by Eq. (87), The additional superscripts o 
and n denote that the radii of curvature are calculated at the re- 
flection boundaries and (2n-l )tt-({i' , respectively. 
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Although the reasoning ertployed to find the distance parameters is 
the same as that used in the py^ceding cases, namely that the total 
field be continuous at the shadow and reflection boundaries, a problem 
arises which was not encountered earlier. For a given aspect of in- 
cidence it is clear that only two of the boundaries associated with the 
three transition functions exist, the other boundary is outside real 
space. Since neither the field or source points are permitted close to 
grazing incidence at = 0 or mr, it is reasonable to set the transition 
function, which is associated with the boundary located outside the 
interval 0 < (ji < ott, equal to one. 

At grazing incidence <()' = tt or (n-l)TT for which or l'^'^ vanish, 
the scalar diffraction coefficients are calculated by the same procedure 
used for the curved screen at grazing incidence (|)* = u. 

In the far zone where s » the principal radii of curvature pp pg 
of the incident and reflected wavefronts at and the radius of curvature 
p of the diffracted wavefront at in the directions of incidence and 
reflection, Eqs. (77), (87a), and (87b) simplify to the form 


(90) 


L = 


P1P2 


sin R 


the appropriate superscripts are omitted here for the sake of riOtational 
simplicity. 

An interesting case occurs if there is a caustic of the incident, 
reflected or diffracted wave on a shadow or reflection boundary. The 
radii of curvature p^, P2 or p associated with such a caustic are negative, 
and L may be either negative or positive. If L is positive, the presence of 
caustics at these boundaries presents no difficulty, except at points near 
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the caustic itself. On the otherhand if L is negative, there is a 
problem because the transition function has two branches each with an 
imaginary argument. We will restrict our attention to the situation where 
all the caustics on the boundary lie between the field point and the edge; 
this may occur in far-zone field calculations for example. 

As pointed out in the last paragraph of Appendix II, if L is negative 
the incident (or reflected) field has one more caustic on the shadow (or 
reflection) boundary than does the diffracted field. This means that the 
phase of the transition function must change by an additional tt/2 as one 
moves from a point outside the transition region to the boundary, so that 
the transition function must have a total phase variation of 3tt/4 instead 
of the tt/4 phase variation shown in Fig. 6, An examination of the two 
branches of the transition function at the boundary and outside the 
transition region reveals that they do not have the proper behavior. 

When a curved strip is illuminated by a plane wave from its concave 
side, there is a caustic of the reflected field on the reflection boundaries. 
In treating the scattering from this strip we have found that an adequate 
function is provided by 

in which F(klL|a) is the ordinary transition function given by Eq. (51). 

(Note that L and a may have superscripts). In spite of the fact that 
the above expression has the proper behavior outside transition regions 
and at shadow or reflection boundaries and also appears to yield good 
numerical results, it lacks theoretical justification. A satisfactory 
derivation of the transition function for L negative is being sought. 
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IV. DISCUSSION 

A dyadic diffraction coefficient has been obtained for an 
electromagnetic wave obliquely incident on a curved edge formed by 
perfectly-conducting curved or plane surfaces. Unlike the edge 
diffraction coefficient of Keller's original theory, this diffraction 
coefficient is valid in the transition regions of the shadow and 
reflection boundaries. Although the diffraction coefficient has 
been given in dyadic form in the earlier chapters, it can also be 
represented in matrix form, so that the high-frequency diffracted 
electric field can be written 



and since the high-frequency diffracted magnetic field 
(92) s X E^, 



in which D^, are given by 

(a) Eq. (89) for the curved edge (general case), 

(b) Eq. (86) for a curved edge in a curved screen, 

(c) Eq. (83) or (84) for a curved or straight edge in a plane screen, 
and p is given by Eq. (32). 
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(d) Eq. (50) for the wedge. 

It is pointed out in Section I I IB that the scalar diffraction 
coefficients in cases (a) and (b) are not valid at aspects of 
incidence and diffraction close to gracing on a convex surface 
forming the edge at the point of diffraction. Work is in progress 
to remove this limitation. Also grazing incidence on a plane 
surface is a special case which requires the introduction of a 
factor of 1/2 when calculating the diffracted field, see the 
discussion of page 27. 

The large parameters are kL or kL , kL in the asymptotic 
approximation; hence when these are small our GTD representation 
of the diffracted field is no longer valid. Thus source or field 
points close to the edge (s or s' small) must be excluded; also 
aspects of incidence close to edge-on incidence (b^ small) must be 
excluded. Edge-on incidence is a separate phenomenon, which has 
been discussed by Ryan and Peters [28] and by Senior [29]. 

Outside of the transition regions where the arguments of the 
transition functions are greater than 10, the expressions for the 
scalar diffraction coefficients all simplify to Eq. (49). Usually 
the field point is only in one transition region at a time, so 
that the calculation of the diffracted field is simplified because 
only one of the transition functions is significantly different from 
unity. 

One would expect the diffraction coefficients for the wedge to 
be more accurate than those for the curved edge, because the 
canonical problems involve wedge diffraction. If the curved edge 
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were used as a canonical problem, one would anticipate the presence 
of additional terms in the asymptotic solution for the diffracted 
field; these terms would depend upon the radius of curvature of 
the edge at the point of diffraction and its derivatives with 
respect to distance along the edge. This is verified by the work 
of Buchal and Keller [30] and Wolfe [31], who treated the diffraction 
of a scalar plane wave normally incident on a plane screen with a 
curved edge. 

In calculating the diffracted field, it is assumed that the in- 
cident field is slowly varying at the point of diffraction, except 
for its phase variation along the incident ray. If the incident field 
is rapidly varying at the point of diffraction, it is usually possible 
to express it as a sum of slowly- varying component fields, so that 
the diffracted field of each component can be calculated in the usual 
way and the total diffracted field obtained by superposition. Al- 
ternatively, in calculating the diffracted field, one could introduce 
higher order terms which depend upon the spatial derivatives of the 
incident field at the point of diffraction. Expressions of this type 
were obtained by Zitron and Karp [32] in their treatment of the 
scattering from cylinders; they are also derived in Reference 11. 

In the text it is pointed out that Eqs. (91) and (93) can not be 
used to calculate the field at a caustic of the diffracted ray. At 
such a caustic it is convenient to use a supplementary solution in 
the form of an integral representation of the field. The equivalent 
sources in this representation are determined from a suitable high- 
frequency approximation, such as geometrical optics or the GTD. In 
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the case of an axial caustic, it is convenient to employ equivalent 
electric and magnetic edge currents introduced by Ryan and Peters [33]; 
the use of these edge currents is also described in Reference [34]. 

In conclusion we note that the geometrical optics field and 
our expression for the edge-diffracted field are both asymptotic 
solutions of Maxwell's equations. The total high-frequency field 
is the sum of these two fields, and away from the edge it is 
everywhere continuous, except at caustics. Our solution reduces 
to known asymptotic solutions for the wedge, and it has been found 
to yield the first two or three terms in the asymptotic expansion 
of the diffracted fields of problems which can be solved differently. 
Furthermore, the numerical results obtained by its application to 
a number of examples are found to be in excellent agreement with 
rigorously-calculated and measured values. Also we have been able 
to show that our solution is consistent with the reciprocity 
principle, see Appendix IV. 
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APPENDIX I 

THE CAUSTIC DISTANCE FOR REFLECTION 


The principal radii of curvature of the reflected wavefront 
p-j, P 2 and the principal directions (axes) of the wavefront will be 
determined in this appendix. The plane of incidence may be different 
from the principal planes of the reflecting surface, so that the 
principal directions of the incident wavefront are quite distinct 
from those of the reflecting surface. 

This problem has been treated both by Fock [35] and by Deschamps; [36] 
however they did not find the principal radii of curvature nor the 
principal directions of the reflected wavefront. Fock used a surface- 
fixed coordinate system to formulate his solution, and he evaluated 
the resulting 3x3 determinant for the divergence factor D(S) of the 
reflected wave after some rather complicated tensor analysis. On the 
otherhand, Deschamps formulated his solution in a ray-fixed coordinate 
system,* and employing elementary matrix theory together with straight- 
forward coordinate transformations, he obtained a 2 x 2 curvature matrix 
for the reflected wave from that of the incident wave. We will find 
p^, P 2 and their principal directions by diagonalizing his curvature 
matri x. 

Let us begin by defining the curvature matrix employed by 
Deschamps. Consider the curved surface 


*The advantages inherent in using ray-fixed coordinate systems in 
treating ray optical problems have already been noted in the text. 
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(A-1) z = f(x,y) 


with the z axis normal to the surface at the point P where x,y,z = 0. 
Thus 


(A-2) 

where 


f = 




f.. = 


_ 8f 


3X * XX 


3^f 


etc. 


are evaluated at x,y = 0. 

In the neighborhood of the origin 


(A-3) 


z = 


ICf 


,, x2 . 2 f^ zy + y2] 


2 2 
LP] P£J 


in which p^, p£ are the principal radii of curvature in the principal 

A ^ 

directions X, Y respectively. With 


x = x X + y y and X = X X + Y Y 


Eq. (A-3) may be written in matrix notation 


(A-4) z = ]-x Q X = jX Qq X, 

in which Q is a 2 x 2 symmetric matrix and is its diagonal form. 
The matrices Q and are referred to as curvature matrices 
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(A-5) Q = 


3 


f f 

XX xy 

f f 

yy yy 


and 

(A-6) Qq = 

The Gaussian curvature of the surface 

(A-7) K = determinant of Q = |Q| = IQ^^I . 

Let a wavefront be incident on a curved surface S at Qj^ as 
shown in Fig. 2A. 

A A 

Ui» U 2 are unit vectors in the principal directions of S at Qj^ 
with principal radii of curvature R-j, R 2 . 

X-|, X 2 are the principal directions of the incident wavefront at 
Q|^ with principal radii of curvature p^, p^. 

Xp x^ are unit vectors perpendicular to the reflected ray; they 
are determined by reflecting the unit vectors Xp X 2 in the plane tangent 
to S at Q|^, i.e., 

(A-8) x!| = X] - 2(n . X])n, 

2 2 2 

see Fig. 2A. As will be seen, Xp x^ are not in the principal directions 
of the reflected wavefront. We now define 
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Fig. lA. A smooth curved surface near the point P. 



Deschamps has shown that the curvature matrix for the reflected wavefront 
(A- 12 ) Q*" = QJ, - 2(0'^)^ cos 0’ 
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w> 



Fig. 2A. Geometry for the analysis of the reflected 
wavefront. The reflecting surface is S. 

Intersection of a principal 

plane of S at Qr with S 

^Intersection of the plane of 

incidence with the plane 
tangent to S at Qp 

Extension of the reflected ray 

below S. 
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in which the superscript -1 denotes the inverse matrix, the superscript 
T denotes the transpose matrix and is the angle of incidence as before. 
A result equivalent to Eq. (12) has been obtained by Fock. 


(A-13) Q*" = 


^12 



^22 


where 


(A-14a) 



2 cos 6^ 



( ©22 ^ ^ ^ 21 ^ 


(A-14b) 


r,r _ 2 cos 0 


® 22®12 

«1 


^ 11®22 


(A-14c) Q^2 = - 


1 

P2 


2 cos 0 


10 


( 012 )^ (011 
”” ' '' “ “rT ~ 


R 


1 


wi th 


(A-Ud) dj, = XJ . U, . 


We have diagonalized Q*" to find its eigenvalues l/p^, l/p^* 
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in which the + sign associated with p!j^ and the - sign with p^. As noted 
in the discussion following Eqs. (23a, b) in the text, this equation has 
the form of an elementary mirror formula, except that the reciprocal of 

the object distance is replaced by the mean curvature of the incident 

/ 

wavefront. 

The incident spherical wavefront is frequently of interest; let 
us simplify Eq. (A-15) for this case. For the spherical wave Xp 

a4 

can be chosen in any way convenient. Let X^ be in the plane of 
incidence; X 2 is then parallel to the plane tangent to the surface of 
reflection at Qj^. Referring to Fig. 2A, 

A«J « A « A ■■ A 

(A-16a) X-| = -cos e sin u + cos 0 cos w U 2 + sin 0 n, 

(A-16b) Xg = cos u) 0^ + sin 0 ) O 2 , 

A A 

where u is the angle between the plane of incidence and the n, U 2 principal 
plane of the reflecting surface. 
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(A-]7a) 


®7T = - cos 0 ^‘ sin 


(t) 


(A- 7 7b) 


®72 = cos 0^ cos 0 , 


(A-7 7c) 


®27 cos 0 ) 


(A-7 7d) 


®22 ■ sin 0) 


Using the above Eqs. 


W-I8) 


cos 0 ^; 


furthermore 


(A-79a) 



(A-79b) 


P] P£ 


1 = 0 . 


•fence, substituting E,s M 771 , 

OS. (A-77), 


(A-I5), 
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(A-20) I- = 


1. 

r 

Pi 

2 


S T 

COS 0 


r 2 2 

sin 6 sin 6 
2 1 
- + — s — 


R 


1 


1 

r 2 - 2 "1 

sin 0 ^ sin 0 
2 1 

2 

-4 

^ cos^e’ 

L "1 . 

R 1 R 2 


with 


2 2 2 ' 1 
(A-21a) sin 0 ^ = cos co + sin^w cos 0 , 


(A-21b) sin^02 = sin^to + cos^oj cos 


0.| is the angle between the direction of the incident ray s^ and D-i, 

^ M A 

whereas 0g is the angle between S and^Ug. Equation (A-20) was obtained 
by KQuyournjian several years earlier using a different method. ^ 

We conclude this section by giving the eigenvectors .of Q^.; 
these yield the principal directions of the reflected wavefront with 
respect to the Xp x^ coordinates. 


(A-22) 


(Q22 ■ r ^ ^1 
or _ L Pi 

^1 


■ ^12 ^ 


JiQ22 “ ^ 

Pi 


- Q?2 


(A-23) Xj = -s’" X Xlf . 
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Fig. 3A. Principal directions of the reflected wavefront with 
respect to the x||’, coordinates. 


It should be noted that the principal directions of the wavefront 
are distinct from the principal directions associated with the reflection 
matrix; as pointed out in the text the latter are parallel and perpendicular 
to the plane of incidence. 
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APPENDIX II 

THE EDGE CAUSTIC DISTANCE 



Fig. 4A, Caustic distances associated with the diffraction 
by a curved edge. 

A wavefront is incident on a curved edge at the point Q^, where 
the radius of curvature is a, the unit tangent vector is e, and the 
unit normal vector to the edge is n^ given by the Serret-Frenet formula 

(A-24) i ^ 

in which the super dot denotes differentiation with respect to the distance 
along the edge and a is taken to be positive. It is apparent from the 
above equation that n^ is directed away from the center of curvature. 
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A second useful formula will be derived next. Let r be the position 


vector from a point 0 to a point on the edge. 
(A-25) e = r 


(A-26) r = r/r 


from which 


(A-27) 



Since 


(A-28) r^ = r • r , 


(A-29) r = r • r , 


where Eq. (A-26) has been used. 

Substituting Eqs. (A-25), (A-26) and (A-29) into E( 
dot multiplying both sides of the resulting equation by 


A- A 


(A- 30) e • r 




According to Keller's law of edge diffraction. 


A A 


(A-31). s' • e = s • e = cos 6- 


and so 


l, (A-27) and 
e, one obtains 
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A AAA A A A A 

(A- 32) s' •e + e*s'=s*e + s*e 


Employing Eqs. (A-24) and (A- 30) in Eq. (A- 32) 


A A 


(A-33) 


1 - (e » s') 


n_ • s 


A A 


_ 1 - (e • s)^ "e • ^ 


where = s' is the radius of curvature of the incident wavefront at 
taken in the plane containing the incident ray and e, 
p is the caustic distance; it is also the radius of curvature 
of the diffracted wavefront at Qg taken in the plane containing 

A 

the diffracted ray and e. 

Using Eq. (A-31), Eq. (A-33) may be written more compactly as 


(A-34) 


A A A 

1 ■ 1 "e • 

0 I , 


which has the same form as the elementary lens equation, where p^ and p 
correspond to the source and image distances, respectively. Thus we 
may write Eq. (A-34) as 


(A- 35) 




where f is interpreted as the focal distance. Eq. (A-34) was first given 
by Kouyoumjian [8], and it has also appeared in a number of ElectroScience 
Laboratory reports [34,37]. Recently Deschamps has rederived Eq. (A-34) 
by a different method [36]. 
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An alternative expression for p has been given by Keller [3]; 
we will derive it here for the sake of conpleteness. 

From Eq. (A-31 ) 

A A A A 

s*e + s*e = -sin 


Again using Eqs. (A-25), (A-30) and (A-31) 
(A- 36) 


2 

sin B^ n„ • d 
0 _ e 

■; P ■ a 


= - sin SoSq . 


(A-37) 

} , 



P ' Sin Bq 


+ 


cos 6 

7 ~ 2 ~ 1* 

a sin Bq 


Here 5 is the angle between n^ and the diffracted ray. 

It has been found more convenient to calculate the edge caustic 
distance p from Eq. (A-34), which does not contain any derivatives. 
Since the edge is a caustic of the diffracted rays, it is clear 


that one of the principal directions of the diffracted wavefront lies 
in the plane of diffraction and is perpendicular to the diffracted ray. 

The other principal direction is perpendicular to the plane of diffraction 
and is thus tangent to the cone of diffracted rays. At a distance s from 
the point on the edge the principal radii of curvature of the dif-' 
fracted wavefront are p + s and s. 

In general, the principal directions of the incident and reflected 
wavefronts associated with rays incident and reflected at (recall that 
these rays lie on the cone of diffracted rays) do not coincide withthe 
principal directions of the diffracted wavefronts associated with rays. 
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diffracted in the directions of incidence and reflection. However in 
the edge-fixed planes of incidence (or reflection) the curvature of the 
diffracted and incident (or reflected) wavefronts is the same. The 
relationship of the diffracted wavefront to the incident wavefront at 
the shadow boundary is apparent from Eq. (A-34); on the shadow boundary 
s = s' and so p = p^. The relationship between the diffracted wavefront 
and the reflected wavefront on the reflection boundary could be deduced 
from Eqs. (A-15) and (A-34), but it may be verified more easily in the 
following way. From Eq. (33) the reflected and diffracted phase 
^ functions are the same along the edge, i.e., 

/(Qg) = ♦‘‘(Qe) ! 


hence 

<P*^(Q^) .+ s = <|;^(Q£) + s 

in the plane containing the tangent to the edge e and the reflected, 
diffracted rays at the reflection boundary. It follows that p = p^, 
the radius of curvature of the reflected wavefront in the plane 
containing the reflected ray and e. The relationships between the radii 
of curvature of the diffracted wavefront and the incident and reflected 
wavefronts at the shadow and reflection boundaries described in this 
paragraph have been noted by Lewis and Boersma [12], 
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APPENDIX III 

THE PLANES OF INCIDENCE AND REFLECTION 


^ ; '.A . . r ' . • • ' ’ • 

In this appendix it will be shown that the angle a between the 
edge-fixed plane of reflection and the ordinary plane of incidence is 
equal to minus the angle between the edge-fixed plane of incidence and 
the ordinary plane of incidence. We will refer to the ordinary plane 
of incidence simply as the plane of incidence. 

The plane of incidence contains the unit vectors s' and n or 

s and n. The edge-fixed plane of incidence contains the unit vectors 

s' and e. The edge-fixed plane of reflection contains the unit vectors 

s*^ and e. The plane of incidence and the edge-fixed plane of incidence 

intersect along the incident ray at Q^. The plane of incidence and 
the edge-fixed plane of reflection intersect along the reflected ray 
at Qg. 

A^ A 

Referring to Figs. 2, 5 and 5A and recalling that s = s' at the 

edge. 


(A-38a) 


gi _ s' X (n X s' ) 

" |s' X (n X s' ) I , 


(A-38b) 5^ 


-s' X (e X s' ) 

A A A 

|s' X (e X s' ) I , 


(A- 38c) 


*r _ s X (n X s ) 

" ' |s»'x (n x^")! , 
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Fig. 5A. Geometry of the planes of incidence; edge-fixed incidence 
and edge-fixed reflection. The incident and reflected 
rays are perpendicular to the page and directed outward. - 


(A-38d) 





Let a be the angle between the plane of incidence and the edge 
fixed plane of reflection, see Figs. 9 and 5A, 


Am a 

COS a = e • gg 


A Am a Am 

-(n • 


Let be the angle between the plane of incidence and the 
edge-fixed plane of incidence. 
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(A- 39) 


1 1 
cos a = e 


g* = (n»s')(e-s') ^ ^ 

f-'r\ A A O AAA O 


Since (n*s‘) = ~(n*s**) from the law of reflection and (e*s‘) = (e-s*^) 
from the law of edge diffraction, it follows from Eqs. (A-38) and (A-39) 
that 


(A-40) ■ cos = cos a; 
hence 


(A-41) ct^ =±a . 

From the definition of angles in Fig. 5A, = -a. 
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APPENDIX. IV . 

RECIPROCITY . : 

The Lorentz reciprocity theorem inposes an important condition. on 
the solution of an electromagnetic problem. This is a necessary con- 
dition and approximate solutions vdiich fail to satisfy it may contain 
errors which are unacceptably large. We will develop the reciprocity 
condition for the edge structure under consideration in this report, 
and- see if the high-frequency solutions given in Chapter HI satisfy’ it. 
Also it will be shown how reciprocity may be used to extend our analysis 
to a special configuration where it would not be valid if applied ■ 
directly. 

Let the sources g and g' radiate in the presence of an edge 
structure as shown in Fig. 6A; the sources are positioned at a finite 
distance from the edge and from each other. The surface S rests on the 

perfectly-conducting boundary, except in the imnediate vicinity of edge, 

/ 

and extends to infinity. The surface S is the s urf a ee^at- infinity; it 

00 

joins with S and together with this surface encloses the region occupied 
by the sources g and g' . The surface S is sharply rounded at the edge 
so as to enclose it tightly, but its radius of curvature does not vanish 
there. 

The sources g and g' consist of electric and magnetic current 
moments which may be distributed in a volume, over a surface or along a 
line, so that the infinitesimal electric current moment 


76 



1 0^ dv, in a volume 

T ds, over a surface 
§ 

Ij, di!,, along a line 


and the infinitesimal magnetic current moment 




Ky dv, in a volume 
(A-42b) ® surface 


Mi, d£, along a line 


in which is the volume density of electric current, is the volume 
density of magnetic current, etc. Starting with Maxwell's curl equations 
in the form 


^ dp. 


(A-43a) V X H - jwe E = 




(A-43b) -V X E + J03y H = 


dp, 


m 




with a little manipulation [38] in which the divergence theorem, is employed, 
one obtains . , 


7 (r • dPg - F' * dp^) - 

source g 

f(F • dF' - TT • dF') = 

j '^e m 

source g' 


where the source g' and its fields are denoted by primes and the source 
g and its fields are the unprimed quantities. 



Fig. 6A. Region pertaining to the discussion on the 
reciprocity principle. 


The fields of g and g' satisfy an edge condition, also n x E, 

A 

n X E' = 0 on the boundary surface; hence there is no contribution to 
the surface integral in Eq. (A-44) from the integration over S. The 
fields of g and g' also satisfy a radiation condition, so the 
integration over in Eq. (A-44) vanishes too. Thus 
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(A-45) 




(P 


dp^ - 


‘^n,) = 

m 


source g 


(r. 

source g' 


dp; - H • dp;) 


which is a statement of reciprocity in integral form. Since the above 
must be true for any distribution of the sources g, g'. 


(A-46) E' • dp^ - H' • dp„ = E • dp’ - H • dp' 

It is well known that the geometrical optics field satisfies the 
reciprocity condition. Since the total field is the sum of the geometrical 
optics field plus the diffracted field it follows that the diffracted 
field also satisfies the reciprocity condition; in other words 


(A-47) E^' • dpg - H^' • dpjj^ = E^ • dp^ - 11^ • dp^^ . 


Now let us see if the 6TD solution for the diffracted field satisfies 
Eq. (A-47). We iVyill first consider the reciprocity relationship for a 
pair of incremental electric current moments d^^, d^^. The electric 
field of d^g incident on the edge is 


(A-48) 



s' X (s' X dpg) 


(A-49) E^ • dpg = a' • D(<t.,(t)' ; 3 ^) • a h(s,s') 
E^ • d^^ = ? • O’ (<{,', • "a' h(s',s) 


(A-50) 



in which 


(A-51a) a = 


jkz 

-rr 


(A-51b) a' = 


jkz. _ 

' = dp' 

J 4tt 


(A-51c) h(s.s')= 


It I ■ P ■ 

J s(p+s) 


(A-51d) h(s',s) = 


1 / p\ , n-jk(s'+s) 

s Js^{p'+s ' ) ® 


(A_52a) 


(A-52b) ^ ^ . I _ 


and 


1 n • (s'-s) 

(A-53) 7 ~ ~ TT 

Pg sin 8 

is unchanged when the source and field points are interchanged. Note 

A A 

that the directions of s' and s are reversed when they are interchanged. 
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Employing Eqs. (A-52) it is seen that 


(A-54) (s'2s)(^) = 

= (ss.)^(Ui, 

MsV)(a^) 

From Eqs, (A-51 ) and (A-54) h(s,s'),has the useful symmetry property 
(A-55) h(s,s' ) = h(s' ,s) . 

To satisfy the reciprocity relationship for electric current sources, 
(A-56) E^' • dpg = E^ • dp' . 

it follows from Eqs. (A-49), (A-50), and (A-55) that it is necessary for 
(A-57) a' • !(<}),((.' ;6 q) • I = F • 'P(<f ' .4;6 q) • F' . 

It is apparent from the form of n in Eqs. (45) and (A-62) that the above 
is true if 

(A-58) Dg(<^,<|)';eQ) = Dg(<^' ,<^;B q). 

h h 

Consequently, to show that our solution satisfies Eq. (A-56), it is 
sufficient to show that Eq. (A-58) holds. Furthermore, this is sufficient 
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to show that the general expression fpr reciprocity given by Eq. (Ar47) 
is also satisfied, because if Eq. (A-58) is true it can be shown in a 
manner analogous to the preceding development that ~ ^ ’ ^^m 

We will examine the scalar diffraction coefficients for wedge , 
diffraction first. Referring to Eq. (50), the quantity 

c°t F[kL + cot^-Sii^i^jF[kL a'(*+*')] 

remains unchanged if (|) and c|)' are interchanged. The quantity 

F[kL cot^ "~|^~'^’J-j F[kL a'(+-+)]. 

also remains unchanged, because the two terms are interchanged if (j> and 
<j)' are interchanged. This can be shown by noting that when <}. and (j>‘ are 
interchanged 

is replaced by - N”, 

N” is replaced by - N^, 

thus a"'^((j>'~(|)) = a~((t>-<t>' ) , 
and a“(<j)'-({i) = a"'^((j)-(|)' ) . 

For the sources dF , d^' 

6 6 

(A-59) L = sin^a^. 
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so that it is invariant with respect to an interchange of points defined 
by ? and . 

From the preceding discussion it is clear that the scalar 
diffraction coefficients for the wedge satisfy Eq. (A-58); hence 
our GTD solution for the diffraction of an electromagnetic wave by 
a perfectly-conducting wedge is consistent with reciprocity both 
within and outside the transition regions. 

Turning next to the scalar diffraction coefficients for the 
curved edge in the curved screen, it is apparent that Eq. (86) is 
unchanged when <j> and (j>' are interchanged. We note that is in- 
variant with respect to an interchange of source and observation 
points, because it is given by Eq. (A-59), Thus if the pair of sources 
are positioned so that they are not in transition regions of each 
others reflection boundaries, F[kL*"a (<{.+((>' )] can be replaced by unity, 
and Eq. (A-58) is satisfied. On the other hand, if the pair of 
sources d^^, dp^ are positioned on each other reflection boundaries, 
the quantities 


(A-60a) 



(A-60b) B = 



= i + :L = i+ 1 + 1 


r 

P2 


(A-60c) C = ^ = 

' ' pS S p s s f 
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are invariant when s and s' are interchanged. Since 
(A-61) L'- = -^sin\. 


it too must be invariant with respect to this interchange, with the result 

that Eq. (A-58) is satisfied. Thus we have been able to show that our 

GTD solution for the diffraction of an electromagnetic wave by a curved 

screen satisfies the reciprocity condition everywhere except in the 

transition region of a reflection boundary, and within this region it 

satifies the reciprocity condition directly on the boundary itself. It 

is therefore reasonable to assume that this solution is consistent with 

reciprocity at all points where it is valid, with the possible ex-^-. 

ception of some small departures from reciprocity at certain points ; 

in the transition region of the reflection boundary.* 

It is pointed out in the text that as (}>' approaches tt and one is 

close to grazing incidence on the curved screen (see Fig. 11), the 

solution for the diffraction coefficients given by Eq. (86) becomes 

invalid, because approaches o. However, we may use reciprocity to 

calculate E^(p) in the following way. Set d^^^ equal to o and replace 

d^g in Eq. (A-47) by a unit electric current moment u directed parallel 

to the component of E^ (p) which is to be calculated; then for the 

incremental sources dF' and dF' at 0 in Fig. 11 

'^e m 

*The solution may satisfy reciprocity at these points, but we have 
not been able to show it. 
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(A-62) u • E^'(p) = 


— A A A A 

in which D((|)' ,<1 );Bq) = -3^3^ Dg(<{)' ,<!>;3 q)- 4><{>' D^(<J>' ,(j>;3Q) 


s) 


and h(s',s) is given by Eq. (A-51d). It is assumed here that kL*^ for 
incidence from P is sufficiently large so that our approximation for 
the diffracted field at 0 is valid. In Eq. (A-62) one would ordinarily 
let 





<t> 




The discussion of reciprocity for the curved edge in the 
otherwise smooth, curved surface follows in the same manner as 
that for the curved edge in the curved screen. One concludes from 
this that Eq. (89) satisfies Eq. (A-58) for all points at which the 
scalar diffraction coefficients are valid, except possibly for some 
points in the transition regions of the shadow boundaries, were small 
deviations from reciprocity may occur. Also if a problem arises in 
calculating the diffracted field near grazing incidence on a curved 
surface (f'^r)* this calculation may be carried out with the aid of 
Eq. (A-62). 
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